A small superconducting grain in the canonical ensemble 
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Abstract 

By means of the Lanczos method we analyze superconducting correlations 
in ultrasmall grains at fixed particle number. We compute the ground state 
properties and the excitation gap of the pairing Hamiltonian as a function 
of the level spacing 5. Both quantities turn out to be parity dependent and 
universal functions of the ratio 5/ A (A is the BCS gap). We then character- 
ize superconductivity in the canonical ensemble from the scaling behavior of 
correlation functions in energy space. 
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What is the size limit for a metal particle to have superconducting properties? An- 
derson [|l] posed this question back in 1959 arguing that when the average level spacing 
S (inversely proportional to volume of the grain) becomes of the order of the BCS gap A 
superconductivity should disappear. A related question is how to characterize " superconduc- 
tivity" in small systems. The transition is washed out, for instance, by thermal fluctuations 
of the order parameter 0. Moreover the hallmarks of Cooper pair condensation like the 
zero resistance and the Meissner effect are absent when the grains are of submicron size. 

In a series of recent experiments Ralph, Black and Tinkham p,|J studied the transport 
through nanometer-scale Al grains. These experiments revealed the existence of a spectro- 
scopic gap larger than the average level spacing which could be driven to zero by applying 
a suitable magnetic field. This was convincingly interpreted as the reminiscence of super- 
conductivity. As the grain size was further reduced ( < 5 nm) no trace of the gap in the 
spectrum was detected. The experimental results were found to be parity dependent, i.e. 
depending on the electron number in the grain being even or odd. 

In the light of these experiments von Delft et al. reconsidered the question posed by 
Anderson. They included the effect of a uniform finite level spacing in a parity dependent 
mean field theory , and found that the breakdown of superconductivity (in the BCS sense) 
occurs at a value of 5/ A which is indeed parity dependent. In grains with an even number 
of electrons superconductivity persists down to smaller grain sizes as compared with the odd 
ones. This parity effect gets enhanced when the effect of level statistics J7| is included. 

Due to the spontaneous breaking of the gauge symmetry, the BCS theory is most trans- 



parently formulated in the grancanonical ensemble |L0| since it is easy to define the order 
parameter as the amplitude to create (or destroy) a Cooper pair in the condensate. In the 
canonical ensemble this quantity vanishes and the characterization of coherence is more dif- 
ficult. Fortunately the use of the grancanonical ensemble is appropriate for many systems 
which are large enough to give small relative fluctuations of the electron number. This con- 
dition is not strictly met in the experiments of Refs. |3],^| where charging effects allow to fix 
the number of electrons in the grain ||. For these systems moreover quantum fluctuations 
of the pairing field may become so large to invalidate the mean field approach. In order 
to characterize the ground-state pair correlations, Matveev and Larkin || then proposed to 
use the parity gap Ap, an experimentally accessible quantity related to the extra ground 
state energy of a system with an unpaired electron. Even the limit of very small grains Ap 
is very sensitive to superconducting fluctuations (in the opposite case it reduces to A). 

Superconductivity in ultrasmall grains requires to study simultaneously the effect of finite 
level spacing and quantum fluctuations at fixed particle number. In this Letter we tackle this 
problem by Lanczos exact diagonalization. We characterize the superconducting correlations 
by studying the parity effect in the ground state (Fig.|I] and Fig.|2|) and the spectroscopic gap 
(Fig.|). Then we address the central question of defining superconductivity at fixed particle 
number by discussing the scaling properties, in energy space, of the pairing model (see Fig.|] 
and Fig.^). 

The BCS pairing Hamiltonian for the small grain is 



H ^ ] Cm cJj )(T C nj( j Ot 5 ^ ] C m _(_C r)lj _C nj _C n) -|_ (1) 
11=1 m,n=l 

<T = ± 

The indices m and n label the single particle energy levels with energy e m and annihilation 
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operator c mjCr . The quantum number a = ± labels time reversed electron states. The 
number of (doubly degenerate) levels is fixed to Q which is twice the Debye frequency u D 
in units of 5. Finally a is the dimensionless BCS coupling constant and 5 is the average 
level spacing (5 ~ 1/JV(0)V, N(0) being the density of states at the Fermi energy and V 
the volume of the grain). Since the Hamiltonian contains only pairing terms, an electron 
in a singly occupied level cannot interact with the other electrons (the unpaired electron is 
frozen). In the following we will use the simplified model with equally spaced single particle 
levels e m = 5 m We will comment later about the effect of level statistics [0. We 

study systems up to Q = 25 at half filling (Q = N) which corresponds to the usual case of 
attractive interaction in a shell | e |< ujd centered at the Fermi energy. 

We first consider the properties of the ground state by measuring the parity gaps f| 

Ap = E 2 n+i — ^ (E 2 n + E2N+2) 

Ap = —E 2 n + - {E2N+1 + E2N-1) 

Here E^ is the ground state energy for a system with N electrons. By increasing the 
level spacing, Ap and Ap behave in a different way. The case ^ < A has been discussed 
in ROB, Ap/A Pd 1 — 5/2 A while Ap has only exponentially small corrections, Ap/A rs 



1 — J 5/ A exp{— 2n5/A}. In the opposite limit, 5 ^> A the behavior of the parity effect 
is dominated by strong superconducting fluctuations H which give logarithmic corrections 
to the non interacting result, i.e. Ap = 5/2ln(5/A) (much larger than the BCS gap at 
the same level spacing). By comparing these two limits Matveev and Larkin concluded 
that a minimum should appear when the level spacing is of the order of A. In Fig. [I] the 
results of the numerical diagonalization for Ap are presented as a function of 5/ A. The two 
asymptotic behaviors discussed above are plotted for comparison. In Fig. |2] we plot Ap/A 
which monotonically increases when the grain size is reduced, as expected ||. Given the 
maximum number of levels we can account for, in the limit 5 <C A we need to consider 
couplings up to a = 0.5. In this regime we use the relation A = uod/2 sinh(l/a) to get the 
comparison of our exact diagonalization results with the asymptotics given in Refs. P,P2"|,|5]|. 



Notice that our data, which refer to systems with different number of electrons (from N = 10 
to N = 25) collapse on a single curve for all values of the ratio 5/ A and not only in the 



limiting cases discussed in Refs. [|||T2|,g. This suggests that Ap/A is a universal function of 
5/ A (a detailed account of this issue will be presented elsewhere [p3]). One consequence of 
that is the systems we consider, although small compared with the superconducting grains 
used in the experiments where iV ~ 10 3 — 10 5 , may capture all the relevant features of the 
model, in particular in the interesting crossover region 5 ~ A. 

Next we study the spectroscopic gap E G between the ground state and the first excited 
many body level [0,|J. In the noninteracting case Eg = 5 whereas in the BCS limit either 
it coincides with 2 A (even-iV grains) or it vanishes, Eg ~ 5 2 /2A (odd-iV grains). The first 
excited state for even-iV grains belongs to the subspace in which two unpaired electrons 
occupy two single-particle levels close to the Fermi energy whereas the relevant subspace for 
odd-iV grains is obtained by moving the unpaired electron to the next single-particle state. 
In Fig. |^ Eq is plotted as function of the grain size. For small grains the effect of pairing 
correlations in the even case is still observable in a rather large range of 5/ A for which odd 
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grains have already reached the asymptotic behavior Eg = 5. The crossover to the "strong 
coupling" regime occurs at values of 5/ A which are different in the odd ad in the even case 
and roughly agree with the mean field critical values determined by von Delft et al. in Ref. 
H, 5/ A ~ 4 (even) and 8/ A ~ 1 (odd). It is interesting to notice that the BCS regime is 
reached at (parity dependent) values of 6/ A <C 1 so there is an intermediate region of values 



5/ A < 1 where BCS theory describes well enough only the ground state properties |T3 
Finally we stress that also for Eg an universal behavior is found. 

To summarize this first part we have shown the full crossover between a " weak coupling" 
regime (very small grains, where fluctuational superconductivity manifests itself via loga- 
rithmic renormalizations) and a "strong coupling" regime (very large grains). This situation 
is reminiscent of the antiferromagnetic Kondo problem. The level spacing S is the low energy 
cutoff which tunes the system through the two regimes. The breakdown of the logarithmic 
renormalization marks the crossover to the superconducting phase. This provides a quanti- 
tative answer to Anderson's question |1J despite of the fact that we have bypassed the very 
problem of defining superconductivity. Following the conventional wisdom up to now we 
meant by superconductivity a regime in which BCS results are qualitatively valid. This is 
not satisfactory in the canonical ensemble since the central quantity, the BCS order param- 
eter, is always zero. In order to characterize superconductivity one has to consider higher 
order correlators. However they will be non zero for generic interaction, even for repulsive 
ones, so it is not straightforward to extract from them a quantity which plays the role of 
the "order parameter". 

Nevertheless a characterization of superconductivity in the canonical ensemble can be 
achieved by studying the scaling of correlations in the energy space: a superconducting 
system displays long range energy correlations. To this end we consider the pseudospin 
representation jTJ]] of the Hamiltonian Eq.fJI]) 



H = J2e n (l-2S z n )-a6 £ S+S~ (2) 

n=l m,n=l 

where 5+ = c^ +c^. and S* = (1/2)(1 — c„ i+ c n>+ — _c n _). Each energy level is represented 
by a site of a fictitious lattice and Eq.(Q) is the Hamiltonian of a one dimensional spin- 1/2 
XY model with long range interaction in a nonuniform transverse field |15[ . The number of 
pairs fixes the total S z and for odd electron numbers one should simply remove the "site" 
occupied by the unpaired electron. In the absence of XY interaction the spins point in the 
z-direction with a domain wall at the Fermi energy separating up and down spin regions. In 
the opposite limit (vanishing transverse field) the system possesses long range order in the 



XY plane fllql - The superconducting properties can be studied by defining the appropriate 

we 



correlation functions in this fictitious lattice (see also Ref. [0). Following Ref. P,|18 
consider the quantity 



^ = J2 U nVn (3) 

n=l 



where v n = (S^S^) 1 ^ 2 and u n = (S'S^) 1 ^ 2 . In the noninteracting case \l> = since both 
u n and v n are step function symmetric around the Fermi energy and hence their product 
is zero. In the limit of very strong interaction the occupation probability for the pairs, 
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as a function of the level position, is roughly uniform. In this case an estimate of the 
energy of this configuration, at half filling, is ~ (1/2 — a) SQ 2 . Note that for a > 1/2 the 
system gains energy, due to pair mixing, from arbitrary high energy levels. So if we enlarge 
the phase space available for coherence (for instance by progressively increasing Q at fixed 
N/Q) a "normal" system will not take advantage from the presence of extra levels whereas a 
"superconducting" does. In other words correlations are short ranged in energy in "normal" 
systems whereas a "superconducting" system displays long range energy correlations. 
The finite size scaling Ansatz for \l/ is 

^ = fl*F((a - a cr )Q 1/u ) (4) 

where a cr is the critical point. In Fig |] and Fig ^ the behavior of \I> is shown for the even and 
the odd case respectively. A scale invariant point is found whose value is parity dependent. 
By collapsing all the data on a single curve, shown for the even and the odd case in the 
inset of Fig ^ and Fig |5| respectively, it is possible to determine the exponent v. We obtain 
a cr = 0.315 ±0.002, rj = 0.94, l/u = 0.26 for the even case and a cr = 0.345 ±0.002, rj = 1.08, 
1/u = 0.35 for the odd case. Both the magnitude of \1/ and the location of the critical point 
(see Fig [| and Fig |) confirm the natural conjecture that superconducting correlations is 
destroyed easier in the odd rather than in the even case. 

We propose that the scaling in energy space of properly defined correlation functions (like 
that defined in Eq.([3|)) can characterize quantitatively superconductivity in the canonical 
ensemble. 

The results presented here for the parity gap and the excitation gap are in a good agree- 
ment with the analytical expressions of Refs. [§,0 . If the approximation of equally spaced 
levels is relaxed, mesoscopic fluctuations are expected to be important in the intermediate 
region 5 ~ A ||. Nevertheless the very existence of the quantum phase transition and the 
scaling in energy space is not questioned since it does depend only on the interplay between 
kinetic energy and pairing interaction. The location of the critical point and the values of 
the exponents will be different. 

The inclusion of the level statistics as well as the role of an applied magnetic 
field [[19,2C,18|] will be the subject of a forthcoming publication [13|. 
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FIG. 1. The parity gap is plotted as a function of (5/A for systems with odd electron number. 
All the sets are evaluated at half filling (f2 = N). In the inset the region of very large level spacing 
is plotted to compare with the asymptotic result of Matveev and Larkin. The parity gap is centered 
around N = 17 (o), N = 19 (O), TV = 21 (□), N = 23 (A). In the limit of very small grains we 
got a quantitative agreement by using Ap ~ 5/2 ln(ct<5/A) with a ~ 1.35 
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FIG. 2. The same as in Fig.l for even number of particles. The parity gap is centered around 
TV = 18 (o), N = 20 (□), N = 22 (O). 
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FIG. 3. The spectroscopic gap Eg is plotted as a function of 5/ A for even and odd particle 
systems at half filling. (O fi = 9, 10; A n = 11, 12; < Q = 13, 14; > = 15, 16; V ^ = 17, 18; 
+ fi = 19, 20; x n = 21, 22; * = 23, 24; o fi = 25). 
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FIG. 4. Data for the pair mixing parameter as a function of a = g/S and f2 = N = even, 
which show clearly a phase transition (N = 8 (o), N = 12 (□), TV = 16 (O), N = 20 (A)). In the 
inset it is shown the data collapse. 
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